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Abstract 

This paper considers competitive Lotka-Volterra population dynamics with jumps. 
The contributions of this paper are as follows, (a) We show stochastic differential 
equation (SDE) with jumps associated with the model has a unique global positive 
solution; (b) We discuss the uniform boundedness of pth moment with p > and 
reveal the sample Lyapunov exponents; (c) Using a variation-of-constants formula for 
a class of SDEs with jumps, we provide explicit solution for 1-dimensional competitive 
Lotka-Volterra population dynamics with jumps, and investigate the sample Lyapunov 
exponent for each component and the extinction of our n-dimensional model. 
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1 Introduction 



The differential equation 

dX{t) 



X{t)[a{t)-b{t)X{t)}, t>0, 



dt 

X(0) =x, 

has been used to model the population growth of a single species whose members usually live 
in proximity, share the same basic requirements, and compete for resources, food, habitat, 
or territory, and is known as the competitive Lotka-Volterra model or logistic equation. The 
competitive Lotka-Volterra model for n interacting species is described by the n-dimensional 
differential equation 



dt 



k(t) -J>;(*)*i(*) 



1,2, ■• • ,n, (1.1) I eq08 



where Xi(t) represents the population size of species i at time t, cii(t) is the rate of growth 
at time t, bij(t) represents the effect of interspecific (if i ^ j) or intraspecific (if i = j) 
interaction at time t, ai(t)/bij(t) is the carrying capacity of the ith species in absence of 
other species at time t. Eq. (11. ip takes the matrix form 

^^ = dizg(X 1 (t),... ,X n (t))[a(t)-B(t)X(t)}, (1.2) gpg 

where 

X = (Xi, ■ • • , X n ) T , a = (ai, ■ • • , a n ) T , B = (bij) nxn . 

There is an extensive literature concerned with the dynamics of Eq. (11 .2\\ and we here only 
mention Gopalsamy [I], Kuang [TJ, Li et al. [9], Takeuchi and Adachi [221 123], Xiao and Li 
[21]. In particular, the books by Gopalsamy [4J, and Kuang [7J are good references in this 
area. 

On the other hand, the deterministic models assume that parameters in the systems are 
all deterministic irrespective environmental fluctuations, which, from the points of biological 
view, has some limitations in mathematical modeling of ecological systems. While, popula- 
tion dynamics in the real world is affected inevitably by environmental noise, see, e.g., Gard 
[2, 3J. Therefore, competitive Lotka-Volterra models in random environments are becoming 
more and more popular. In general, there are two ways considered in the literature to model 
the influence of environmental fluctuations in population dynamics. One is to consider the 
random perturbations of interspecific or intraspecific interactions by white noise. Recently, 
Mao et al. [13] investigate stochastic n-dimensional Lotka-Volterra system 

dX(t) = diag(Xx(t), • • • , X n {t)) [(a + BX{t))dt + aX(t)dW(t)\ , (1.3) [jjjg 

where W is a one-dimensional standard Brownian motion, and reveal that the environmental 
noise can suppress a potential population explosion (see, e.g., [Ml [15] among others in this 



2 



connection). Another is to consider the stochastic perturbation of growth rate a(t) by the 
white noise with 

a(t) -> a{t) + a{t)W{t) } 

where W(t) is a white noise, namely, W(t) is a Brownian motion defined on a complete 
probability space (Q, J 7 , P) with a filtration {J 7 }^ satisfying the usual conditions (i.e., it 
is right continuous and increasing while J-q contains all P-null sets). As a result, Eq. (jl.2p 
becomes a competitive Lotka-Volterra model in random environments 

dX(t) = diagpTi(t), • • • , X n {t)) [(a(t) - B(t)X(t))dt + a(t)dW(t)} . (1.4) f^T 



There is also extensive literature concerning all kinds of properties of model f jl.4j) . see, e.g., 
Hu and Wang [5J, Jiang and Shi [6], Liu and Wang [II] . Zhu and Yin [251 126], and the 
references therein. 

Furthermore, the population may suffer sudden environmental shocks, e.g., earthquakes, 
hurricanes, epidemics, etc. However, stochastic Lotka-Volterra model (11.41) cannot explain 
such phenomena. To explain these phenomena, introducing a jump process into underlying 
population dynamics provides a feasible and more realistic model. In this paper, we develop 
Lotka-Volterra model with jumps 

dX(t) = diagpfi(r), ■ • ■ ,X n (r)) \(a(t) - B(t)X(t))dt 

(1.5) |eq42 

+ a(t)dW(t) + I 7(t,u)N(dt,du) 



Here X, a, B are defined as in Eq. (jl.2p . 

a = (a 1: ■ ■ ■ , a n ) T , 7 = (71, • • • , 7 n ) T , 

W is a real-valued standard Brownian motion, N is a Poisson counting measure with char- 
acteristic measure A on a measurable subset Y of [0, 00) with A(Y) < 00, N(dt,du) : = 
N(dt,du) — \{du)dt. Throughout the paper, we assume that W and iV are independent. 

As we know, for example, bees colonies in a field J2D] • In particular, they compete for 
food strongly with the colonies located near to them. Similar phenomena abound in the 
nature, see, e.g., [21]. Hence it is reasonable to assume that the self- regulating competitions 
within the same species are strictly positive, e.g., [251 [26]. Therefore we also assume 



(A) For any t > and i,j — 1, 2, ■ • • ,n with i 7^ j, a«(t) > 0, bu(t) > 0, %(t) > 0, o"j(t) and 
ji(t,u) are bounded functions, ba := un7 g R + bait) > and ji(t,u) > —l,u G Y. 

In reference to the existing results in the literature, our contributions are as follows: 

• We use jump diffusion to model the evolutions of population dynamics; 

• We demonstrate that if the population dynamics with jumps is self- regulating or com- 
petitive, then the population will not explode in a finite time almost surely; 
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• We discuss the uniform boundedness of p-th moment for any p > and reveal the 
sample Lyapunov exponents; 

• We obtain the explicit expression of 1-dimensional competitive Lotka-Volterra model 
with jumps, the uniqueness of invariant measure, and further reveal precisely the sam- 
ple Lyapunov exponents for each component and investigate its extinction. 



2 Global Positive Solutions 



As the ith state X{(t) of Eq. (11.51) denotes the size of the ith species in the system, 
it should be nonnegative. Moreover, in order to guarantee SDEs to have a unique global 
(i.e., no explosion in a finite time) solution for any given initial data, the coefficients of the 
equation are generally required to satisfy the linear growth and local Lipschitz conditions, 
e.g., [T5J. However, the drift coefficient of Eq. (jl.5p does not satisfy the linear growth 
condition, though it is locally Lipschitz continuous, so the solution of Eq. (11.51) may explode 
in a finite time. It is therefore requisite to provide some conditions under which the solution 
of Eq. (II. 5p is not only positive but will also not explode to infinite in any finite time. 

Throughout this paper, K denotes a generic constant whose values may vary for its 
different appearances. For a bounded function v defined on R + , set 

v := inf u(t) and z> := sup v{t). 



solution 



For convenience of reference, we recall some fundamental inequalities stated as a lemma. 
Lemma 2.1. 



x r < 1 + r(x — 1) 



x > 0, 1 > r > 0, 



n 



(l-f)A0|„|p 



X\ P < J2 X i ^ 



n 



(1-§)V0|_|j> 



i=i 



where 



{x G 



Xi > 0, 1 < i < n}, and 
In x < x — 1, 



x\ p , Vp > 0, x G 



x > 0. 



(2.1) leqlOO 

(2.2) [eqlOl 

(2.3) Ieql02 



Theorem 2.1. Under assumption (A), for any initial condition X(0) = Xq 6 
has a unique global solution X(t) e for any t > almost surely. 



1, Eq. (O 



Proof. Since the drift coefficient does not fulfil the linear growth condition, the general 
theorems of existence and uniqueness cannot be implemented to this equation. However, it 
is locally Lipschitz continuous, therefore for any given initial condition X(0) G R™ there is 
a unique local solution X(t) for t G [0,r e ), where r e is the explosion time. By Eq. (jl.5p the 
ith component X^t) of X(t) admits the form for i — 1, • • • , n 



dXi(t) = Xi{r 



di(t) -^bijfyXjtyndt + (Ti(t)dW(t) + I ji(t,u)N(dt,du) 
i=i 
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Noting that for any t G [0, r e ) 

^ (s) _^^ (s)Xj ( s )--a, 2 ( S ) 

3=1 

+ y (ln(l + Ji(s,u)) - ji(s,u))X(du)^ds 

+ [ cTi(s)dW{s)+ [ [ ln(l +ji(s,u))N(ds,du)\, 



together with Xj(0) > 0, we can conclude Xi(t) > for any t G [0, r e ). Now consider the 
following two auxiliary SDEs with jumps 



(2.4) eql03 



dYi(t) =Y i (r)^a i (t)-b ii (t)Y i (t))dt + a i (t)dW(t) + J 7i (t,u)N(dt, du) 
F,(0) =X i (0), 

and 

dZiit) = Zi(t-)[(ai{t) -Y^hA^Yjit) -b u (t)Z i (t)yt + a i (t)dW(t) + J <yi(t,u)N{dt,du) 
Zi(0) =Xi{0). 

(2.5) I eqll2 

Due to 1 + 7i(t, m) > by (A), it follows that for any X2 > x\ 

(1 + ji(t, u))x 2 > (1 + 7»(t, u))xi. 
Then by the comparison theorem [TTJ Theorem 3.1] we can conclude that 



Zi(t) < Xi(t) < Yi(t),t G [0, T e ). (2.6) eqlll 



By Lemma [4.21 below, for 1^(0)(= Xi(0)) > 0, we know that Yi(t) will not be expolded in 
any finite time. Moreover, similar to that of Lemma [4.21 below for Zj(0)(= Xi(0)) > 0, we 
can show 

F(Zi(t) > on t G [0,r e )) = 1. 

Hence r e = oo and Xj(t) > almost surely for any t G [0,oo). The proof is therefore 
complete. □ 



3 Boundedness, Tightness, and Lyapunov-type Expo- 
nent 

In the previous section, we see that Eq. (jl.5p has a unique global solution X(t) G for 
any t > almost surely. In this part we shall show for any p > the solution X(t) of Eq. 
(II. 5p admits uniformly finite p-th moment, and discuss the long-term behaviors. 
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te moment Theorem 3.1. Let assumption (A) hold. 



(1) For any p G [0,1,] there is a constant K such that 

sup E\X(t)\ p < K. 



(3.1) eq9501 



(2) Assume further that there exists a constant K(p) > such that for some p > l,t > 
0, i = 1, - - - , n 



/ | 7 ,(t,M)| p A(rfn) < K(p). 
Jy 

Then there exists a constant A(p) > such that 

sup E\X(t)\ p < K{p). 



(3.2) eq90 



(3.3) eq95 



Proof. We shall prove (13. 3p firstly. Define a Lyapunov function for p > 1 

■n 

V(x) := J2 xP i' x e M +- 

i=l 

Applying the Ito formula, we obtain 

E(e t V(X(t))) = V(x )+E f e s [V(X(s)) + CV(X(s), s)]ds, 

Jo 



(3.4) eq43 



where, for x G Wi and t > 0, 

n 

£V(x,t) :=p^ 



i=i 
n 



X; 



+ E / [(l + 7<(*,t*))P-l-p7^,«)]A(d«)af. 

i=i 

By assumption (A) and ( 13. 2ft . we can deduce that there exists constant K > such that 

n 



(3.5) eq50 



i=l 



- o*P' + ( + p(p^M) x p 



n „ 

+ E / [(i+^u)r-i- Plt (t,u)}\(du)xi 



< A. 



Hence 



E(eV(X(t))) < 7(ar ) + / Ae s rfs = V(ar ) + A(e* - 1), 



which yields the desired assertion (13. 3p by the inequality 
For any p G [0, 1], according to the inequality ( 12, ip . 



exm 



[(1 + 7i (t, u)) p - 1 - pjiit, u)] X(du) < 0. 

Jy 

Consequently 

n 

V(x) + CV(x,t) < [-pbu(tK +1 + (l+pa^t)) 4] , 

i=l 

which has upper bound by (A). Then (13.11) holds with p G [0, 1] under (A). □ 

Corollary 3.1. Under assumption (A), there exists an invariant probability measure for 
the solution X(t) of Eq. (Oj) . 

Proof. Let ¥(t,x,A) be the transition probability measure of X(t,x), starting from x at 
time 0. Denote 

1 f T 

IPr(A):=- J ¥{t,x,A)dt 

and B r := {x G M™ : \x\ < r} for r > 0. In the light of Chebyshev's inequality and Theorem 
133] with p G (0, 1), 

fi T (B c r ) = ±J T F(t,x,B')dt <-Lj\\X(t,x)\vdt < £ 

and we have, for any e > 0, iit{B t ) > 1 — e whenever r is large enough. Hence {ht, T > 0} is 
tight. By Krylov-Bogoliubov's theorem, e.g., [T9l Corollary3.1.2, p22], the conclusion follows 
immediately. □ 

Definition 3.1. The solution X(t) of Eq. (I1.5P is called stochastically bounded, if for any 
e G (0, 1), there is a constant H := H[e) such that for any Xq G K" 

limsupP{|X(t)| < H} > 1 - e. 

t— >oo 



undedness 



As an application of Theorem I3.1[ together with the Chebyshev inequality, we can also 
establish the following corollary. 

Corollary 3.2. Under assumption (A), the solution X(t) of Eq. (II. 5p is stochastically 
bounded. 

For later applications, let us cite a strong law of large numbers for local martingales, e.g., 
Lipster [10], as the following lemma. 



7 



8 numbers 



Lemma 3.1. Let M(t),t > 0, be a local martingale vanishing at time and define 

f* d(M)(s) 

where (M)(t) := (M,M)(t) is Meyer's angle bracket process. Then 

M(t) 



lim 

t— >oo t 



a.s. provided that lim pm(^) < oo a.s. 

t— yoo 



Remark 3.1. Let 



^j oc := | $(t,z) predictable J j |*(s, z)\ 2 X(du)ds < oo 



and for * G ^| 



M(t) 



./Y 



Then, by, e.g., Kunita [HJ Proposition 2.4] 



(M)(t) 



|^(s,^)| 2 A(dw)ds and [M](t) = / / |tf(s, ^)| 2 iV(ds, du), 
'o jy Jo Jy 

where [M](£) := [M,M](t), square bracket process (or quadratic variation process) of M(t). 

Theorem 3.2. Let assumption (A) hold. Assume further that for some constant 5 > — 1 
and any t > 

Ji(t, u) > 8, u G Y, i — 1, • • • , n, (3.6) eq78 

and there exists constant K > such that 

-t r 

\-f(s,u)\ 2 X(du)ds < Kt 



JY 



(3.7) eq99 



Then the solution X(t),t > 0, of Eq. ( 11. 5ft has the property 

1 



lim sup 

t-yoo t 



min b u ft 

H\X(t)\)+ 1J ^- / \X(*)\ds 
n Jo 



< max di, a.s. 

l<j<n 



(3.8) eql6 



Proof. For any x G M™ , let V(x) = x,- n by Ito's formula 



i=i 



ln(V(X(t))) < ln(V(x )) + £ (x T (s)(a(s) - B(s)X(s)/V(X(s)) 

-(X T (s)a(s)) 2 /(2V 2 (X(s))))ds 



+ / X T (s)a(s)/V(X(s))dW(s) + / / ln(l + H(X( S -), s, u))N{ds, du), 
Jo Jo Jy 
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where 

H(x,t,u)= ^2ji(t,u)Xj\ jv{x). 

Here we used the fact that 1 + H > and the inequality (12. 3p . Note from the inequality 
(12.21) and assumption (A) that 

X T (s)(a(s) - B(s)X(s))/V(X(s)) - (X T (s)a(s)) 2 /(2V 2 (X(s)) 

n n n 

E OitfXiis) E Xi(s) E bij(s)Xj(s) 
^ 1=1 1=1 J=l 



— n n 

E^( S ) E^(^) 

i=l i=l 

< max a, =^ Am . 

l<i<n Jn 



t ft 

T 



Let 

M(t):= I X T (s)a(s)/V{X(s))dW{s) and M(t) := / / ln(l + F(X(s~), s, w))iV(ds, dw). 

Compute by the boundedness of a that 

(M)(t) = [\x T (s)a(s)) 2 /V 2 (X(s))ds < f \a(s)\ 2 ds < Kt. 
Jo Jo 

On the other hand, by assumption (13. 6p and the definition of H, for x G R+ we obtain 

H(x,t,u) > 5 

and, in addition to (12.31) . for —1 < 5 < 

| ln(l + H(x,t,u))\ < | ln(l + H(x,y,u))I {5 < H (x,t,u)<o}\ + I hi(l + H(x,y,u))I{ <H(x,t,u)}\ 
<-\n(l + 6) + \H(x,t,u)\. 

This, together with (13.7j) . gives that 

* /• 

2 - 



(M)(t) = / / (ln(l + #(A>),s,u))) 2 A(du)ds 



JY 



< 2(-ln(l + 5)) 2 A(Y)t + 2 f [ H 2 (X(s), s,u))X(du)ds 

Jo Jy 



<2(-ln(l + 5)) 2 A(Y)t + 2 f [ 

Jo Jy 

< (2(-ln(l + 5)) 2 A(Y) + ir)t. 



\j(t, u)\ 2 X(du)ds 

o Jy 



Then the strong law of large numbers, Lemma 13.11 yields 

-M(t) ->■ a.s. and -M(t) ->■ as t -> oo, 
and the conclusion follows. □ 



4 Variation-of- Constants Formula and the Sample Lya- 
punov Exponents 



In this part we further discuss the long-term behaviors of model ( 11. 5p . To begin, we 
obtain the following variation-of-constant formula for 1-dimensional diffusion with jumps, 
which is interesting in its own right. 



t formula 



4.1 Variation-of-Constants Formula 

Lemma 4.1. Let F,G,f,g : R + — > R and H, h : R + x Y — > R be Borel-measurable and 
bounded functions with property H > — 1, and Y(t) satisfy 



dY(t) = [F(t)Y(t) + f(t)]dt + [G(t)Y(t) + g(t)]dW(t) 

+ [[Y(r)H(t,u)+h{t,u))N(dt,du), 
Jy 

Y(0) = Y . 

Then the solution can be explicitly expressed as: 

Y(t) = m(Yo+ /'•-'w [(/(.) -GW«w- f H !- S ;t {s :t' 

V J LV J Y l + H(s,u) 

h(s, u) 



\{du) )ds 



+ g(s)dW(s) 



Y 1 + H(s, u 



-N(ds, du) 



where 



:= exp [ / (f{s) - \g 2 {s) + / [ln(l + H{s, u)) - H( 

+ I G{s)dW{s)+ [ [ ln(l + H{s,u))N{ds,du) 
Jo Jo Jy 



s, u)]X(du) )ds 



is the fundamental solution of corresponding homogeneous linear equation 

dZ(t) =F(t)Z(t)dt + G(t)Z(t)dW(t) + Z(t-) [ H{t,u)N{dt,du). 



Proof. Noting that 

$(t)=exp\ f (F{s) -^G 2 (s) + [[ln(l + H(s,u))-H(s,u)}X(du))ds 
'-Jo ^ 2 Jy ' 

+ / G(s)dW(s)+ [ [ + H(s,u))N(ds,du) 
Jo Jo Jy 

is the fundamental solution to Eq. (14. 2p . we then have 

d$(t) = F(t)$(t)dt + G(t)$(t)dW(t) + $(r) / H(t,u)N(dt,du). 



(4.1) fi^3" 



(4.2) [^8 



(4.3) [e^4 
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By [HE Theorem 1.19, plO], Eq. ( 14.1 p has a unique solution Y(t),t > 0. We assume that 



Y(t) = m (y(o) + J 



f(s)ds + g{s)dW{s) + / h(s,u)N(ds, du) 



where f,g, and h are functions to be determined. Let 



Y(t) = Y(0) + / $- x (s) 



/(s)cfs + ^(s)dW(s) + / h(s,u)N(ds,du) 



which means 



dY(t) = ^-\t) 



f(t)dt + g(t)dW{t) + / h(t,u)N(dt,du) 



(4.4) eq5 



Observing that $ and Y are real-valued Levy type stochastic integrals, by Ito's product 
formula, e.g., [H Theorem 4.4.13, p231], we can deduce that 



dY(t) = $(r)dY(t) + Y(t-)d$(t) + d[$,Y](t), (4.5) [eq6 

where [$, F] is the cross quadratic variation of processes $ and Y, and by (4.14) in |TJ p230] 

d[$,Y](t) = G(t)g(t)dt+ [ H(t,u)h{t,u)N(dt,du). 



(4.6) eq7 



Putting (tO)| . gU), and flU} into (l43j) . we deduce that 



dy(t) 



/(t)dt + g{t)dW{t) + / h(t,u)N(dt,du) 



+ F(t)Y(t)dt + G(t)Y(t)dW(t) + Y(r) / H(t,u)N(dt,du) 



+ G(t)g(t)dt + / H(t,u)h{t,u)N(dt,du) 



f(t) + F(t)Y(t) + G(t)g(t)+ / H{t,u)h{t,u)X{du 



dt + [g(t) + G(t)Y(t)]dW(t) 



Setting 



and 



[h(t, u) + Y(r)H(t, u) + #(f, u)/i(f, u)] iV(df, du). 

f(t) + G(t)g(t) + [ H(t,u)h(t,u)\(du) = f(t) 
Jy 

g(t) = g{t) and h(t, u) + H(t, u)h(t, u) = h(t, u), 



hence we derive that 

/(f) = f(t) - G(t)g(t) - [ ^ff^ A((fa),g(t) = <?(f) and 



/i(t, u) 



1 + H(t,u 
and the required expression follows. 



1 + H(t,u) 



□ 
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4.2 One Dimensional Competitive Model 



solution 



In what follows, we shall study some properties of processes Yi(t) defined by (12 .4p . which 
is actually one dimensional competitive model. 

Lemma 4.2. Under assumption (A), Eq. ( 12 .4p admits a unique positive solution Yi(t),t > 
0, which admits the explicit formula 



Yi(t) 



Toy + J^iWM*)^' 



(4.7) eq9 



Xi(p) 



where 



ds 



:= exp ( j di(s) - ^K s ) + J (M 1 + Ti( s > u )) ~ 7*(*> u))X(du) 

+ [ a l {s)dW(s)+ [ [ ln(l + ^{s,u))N{ds,du)). 
Jo Jo Jy ' 

Proof. It is easy to see that is integrable in any finite interval, hence Yi(t) will never 
reach 0. Letting Yi(t) : = and applying the Ito formula we have 

= -^7^W)[«t) -bu^Yiit^dt + a^dWit)] + l^7^o-!(t)Y?(t)dt 



+ 
+ 



(1 + 7< («,«)W) W) Y?(t) 
1 1 



\{du)dt 



(l + jt(t,u))Yi(t-) Yi(t-) 



N(dt,du), 



that is, 

rff(t) = ?(r)[(V J 2 (t)-a i (t) + 

+ 



- l + 7 i (t,u))A(dM))dt-c7 i (t)dW(t) 



1 



1 + 7i(t,w) 



1 + %(t,u) 
-l)N(dt,du) +b u (t)dt. 

(4.8) [jjg 

By Lemma I4.1[ Eq. ( 14. 8ft has an explicit solution and the conclusion ( 14. 7p follows. □ 

Definition 4.1. The solution of Eq. (12 ,4p is said to be stochastically permanent if for any 
e G (0, 1) there exit positive constants Hi := Hi(e) and H 2 '■— H 2 (e) such that 

liminf ¥{Yi(t) < H x } > 1 - e and liminf P{I$(t) > H 2 } > 1 - e. 



i— >-oo 



t— >oo 



permanent Theorem 4.1. Let assumption (A) hold. Assume further that there exists constant c x > 
such that, for any £ > and i = 1, - ■ ■ ,n, 

Oi(t) -o-iit) - 



-X(du) > ci, 
1 + 7i(t,tt) 

then the solution Yi(t),t > of Eq. ( 12.4j) is stochastically permanent. 



(4.9) eq03 
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Proof. The first part of the proof follows by the Chebyshev inequality and Corollary 13.21 
Observe that (14.71) can be rewritten in the form 



ds 



a i( s ) ~ 2 Cr ^ S ' ) + ' ^ ln '- 1 + ' y i( s > u )) ~ li(s,u))X(du) 



- I ai(s)dW(s)- [ [ ln(l + 7 i (s,u))JV(ds,du)) 







JO J Y 
t , x 



a i r ) - §°"» 2 ( r ) + / ( ln ( 1 + 7*0", «)) ~ li(r,u))\(du) 



+ / 6u(s)exp 
Jo 

-f «r m -r /Ml + *<r,»)W*,*,))* 




dr 



(4.10) f^q02" 



By, e.g., [U Corollary 5.2.2, p253], we notice that 



s,u) 



1 + ln(l + ji(s,u)))X(du)ds 



o-i{s)dW{s) 




ln(l + ji(s, u))N(ds, du) 



is a local martingale. Hence letting Mi(t) := and taking expectations on both sides of 
fl4~T0l leads to 



Xi(0) 



exp 



+ / bu(s)exp 



ai{s) - a^s) 
t 



1 +7i(s,M 



-A(dit) 



a; r - <x, r 



1 + 7i(r,w) 



A(dit) 



drds, 



which, combining (14. 9p . yields 



EM,(f) < Y7m e " C1< + fUs)e- c ^- s) ds <-+ - - 

X t {0) Jo ci \^(°) c i 



-cit 



Hence there exists a constant > such that 

EMi(t) < K. 



(4.11) I eq05 



(4.12) [^M 



Furthermore, for any e > and constant H 2 (e) > 0, thanks to the Chebyshev inequality and 
(ED 

¥{Yi(t) >H 2 } = F {Mi(t) < 1/H 2 } = 1 - P {M t (t) > l/# 2 } > 1 - H^M^t) > 1 - e 
whenever H 2 = e/K, as required. 

□ 
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Theorem 4.2. Let the conditions of Theorem 14.11 hold. Then Eq. (12.41) has the property 

lim E\Yi(t,x) — Yi(t,y)\^ = uniformly in (x,y) 6lxK, (4.13) I eq06 

where K is any compact subset of (0, oo). 

Proof. By the Holder inequality 

E\Y i (t,x)-Y i (t,y)\^E[Y i (t,x)Y i (t,y) ' 1 



< (E(Yi(t,x)Yi(t,y)))2 E 



Yi(t,y) Yi(t,x) 
1 



Yi(t,y) Yi(t,x) 



To show the desired assertion it is sufficient to estimate the two terms on the right-hand side 
of the last step. By virtue of the Ito formula, 

d(Xi(t, x)Yi(t, y)) = Yi(r, x)dY i (t, y) + Y(r ,y)dY i (t, x) + d[Y(t, x), Y(t, y)} 

{a i (t)-b ii (t)Y i (t,y))dt + a i (t)dW(t) + I ji{t,u)N(dt,du) 



f Mt) - bii{t)Yi(t, x))dt + <ji{t)dW{t) + / 7i(t, u)N{dt, du) 



Y i (r,x)Y i (r,y) 

+ Y i (r,x)Y i (r,y) 

+ a*(t)Y i (t,x)Y i (t,y)dt+ I ^(t,u)Y i (r,x)Y i (r,y)N(dt,du) 

Jy 

(2ai(t) + o?(t))Yi(t, x)Yi{t, y)dt - 6«(*)*i(*, x)Y(t, y)(Yi(t, x) + Y(t, y))dt 

+ 2a i (t)Y i (t,x)Y i (t,y)dW(t) + 2 [ j i (t,u)Y i (t-,x)Y i (r,y)N(dt,du) 

Jy 

+ / ^(t,u)Y l (r,x)Y i (t-,y)N(dt,du). 



Thus, in view of Jensen's inequality and the familiar inequality a + b > 2y ah for any a, b > 0, 
we deduce that 

E(Y i (t,x)Y i (t,y))<xy + [ 8 i (s)E(Y i (s,x)Y i (s,y))ds 



E / b ll (s)(Y l (s,x)Y l (s,y)(Y l (s,x) + Y i (s,y)))ds 



<xy+ 5 i {s)E(Y i {s,x)Y t {s,y))ds- / b ii (s)(E(Y i (s,x)Y i (s,y)))lds, 
Jo Jo 

where 5i(t) := 2di(£) + of (t) + J Y jf(t,u)X(du). By the comparison theorem, 



E(Y i (t,x)Y i (t,y))< (l/^ye-yo 5 ^ ds + l f h^e^^^^ds 



(4.14) eq07 



< [b/$i + (l/Jxy-b u /6i)e ^ 
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On the other hand, thanks to (14. 7p we have 



Yi(t,x) Yi(t,y) 



1 1 

j exp 

x y 

t 

ai(s)dW(s 
o 



ai(s) - -a { (s) + j (ln(l +7»(s, u)) - ii(s,u))\(du) 



ln(l + %(s,u))N(ds,du) ). 



ds 



JY 



In the same way as (14. lip was done, it follows from (14. 9p that 

1 1 



E 





1 


1 


< 








X 


y 



C2t 



Yfax) Yi(t,y) 

Thus (Q3l) follows by combining ( Olj) and (OBI . 

If Oj, cjj, 7j are time-independent, Eq. (12 ,4p reduces to 



(at - biiYi(t))dt + <7idW(t) + / ii(u)N(dt, du) 



(4.15) |eq!9 
□ 

(4.16) f^q2T 



with original value x > 0. Let p(t, x, dy) denote the transition probability of solution process 
Y{(t, x) and F(t,x,A) denote the probability of event {Yi(t,x) G A}, where A is a Borel 
measurable subset of (0, oo). It is similar to that of Corollary 13. 1[ under the conditions 
of Theorem 14.11 there exists an invariant measure for Y{(t,x). Moreover by the standard 
procedure [T5], p213-216], we know that Theorem 14.21 implies the uniqueness of invariant 
measure. That is: 



tribution 



Theorem 4.3. Under the conditions of Theorem 14.11 and 14.21 the solution Yi(t,x) of Eq. 
(14.161) has a unique invariant measure. 

We further need the following exponential martingale inequality with jumps, e.g., [TJ 
Theorem 5.2.9, p291]. 



artingale Lemma 4.3. Assume that g : [0, oo) — » R and h : [0, oo) x Y — > R are both predictable 
y^-adapted processes such that for any T > 



/ \g(t)\ 2 dt < oo a.s. and / / \h(t, u)\ 2 X(du)dt < 
Jo Jo Jy 

Then for any constants a, (3 > 



oo a.s. 



P<{ sup 

0<t<T l Jq 



g{s)dW{s) 



a 



\g{s)\ 2 ds + 



h(s, u)N(ds, du) 



JY 



1 

a 



JY 



,«h(.,«) _ j _ a h(s,u)]X(du)ds > /?} < e" Q/3 . 
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property Lemma 4.4. Let assumption (A) hold. Assume further that for any t > and i — 1, 



n 




sup 

t>o Jo Jy 



e s [7j(s, u) — ln(l + ji(s, u))]X(du)ds < oo. 



(4.17) eq32 



Then 



limsup — — - < 1, a.s. for each i — 1, • • • ,n, 
t-^oo m t 



Proof. For any t > and z = 1, • • • , n, applying the Ito formula 



e t lnY l {t) = \nX t (0) + / e s kiY^s) + a^s) - bM^s) - -af(s) 



ds 



+ / [ln(l + 7i(s,u)) - %(s,u)]X(du) 



+ / e s ai(s)dW(s) 




e s ln(l + 7i(s, u))N(ds, du). 



</¥ 



Note that, for c.x > 0, lnx — cx attains its maximum value —1 — In c at x = -. Thus it 
follows from the inequality (12. 3 1) that 



e* In y;(t) < lnXi(O) + / e 



1 - lnfeii(s) + ai(s) - ^K s ) 



ds 



e s o-i(s)dW(s) + 
In the light of Lemma [4.31 for any a, (3, T > 
p{ sup 




(4.18) eql09 



e s ln(l + 7i(s, u))N(ds, du). 



o J Y 



o<t<r 
1 
a 



a 



eV^jdlf (s) - ^ / e 2 V 2 (s)ds + 




JY 




./Y 



3 ae s ln(l+7i(s,u)) 



e s ln(l + 7i(s, u))N(ds, du) 
1 - «e s ln(l + 7i(s,M)) X(du)ds > /?} < e" 



-a/9 



Choose T = &7, a = ee fe7 , and /? 



6»e fc T lnfc 



where fc e N, < e < 1, 7 > 0, and 6 > 1 in 



the above equation. Since ^ A; e < 00, we can deduce from the Borel-Cantalli Lemma that 



k=l 



there exists an C £1 with P(Oi) = 1 such that for any e G fi ; an integer fcj = ki(oj, e) can 
be found such that 



e s o-i(s)dW(s) + 

'0 JY 

delink ee- k ~< rt 
< + 




+ 



e 2 
1 '* 



ee 



-&7 




JY 



e s ln(l + 7i(s, u))N(ds, du) 
e 2s a*(s)ds 

(1 + Ti (s, u)) eeS ^ 7 - 1 - ee s ~ kl ln(l + 7i (s, u))l \{du)ds 







16 



whenever k > ki, < t < kj. Next, note from the inequality (12.1 p that, for any wgQ; and 
0<e<l,0<t<Jfe7 with k > h, 



1 



ee 



JY 
t 



(1 + ji{s, u)) eeS ~ kl - 1 - ee 3 -^ ln(l + 7i (s, u)) A(dw)ds 



s—k'y - 



< / / e s -'( 7i (s,M)-ln(l + 7 i (s,u)))A(^)ds. 
Jo if 

Thus, for a; G fij and (A; — 1)7 < t < k^ with k > k{ + 1, we have 



lnr,(t) < lnX t (0) | 



6e k ^ In fc 



hit ~ eHnt ee( fc - 1 ^ln((A; - 1)7) 
1 '* 



f in7./ ( 



- 1 - lnfefs) + ads) - -(1 - ee s - k ^)af(s) 



ds 



+ — I I e s - t [ li (s,u)-\n(l + li (s,u))]X(du)ds. 
hit Jo Jy 

Letting k "f oo,together with assumption (A) and (14.17[) . leads to 

.. \nYi(t) 9ei 
hmsup — < , 

t^oo mt e 

and the conclusion follows by setting 7 4. 0, e f 1, and 6^1. 

hit 

Noting the limit lim = 0, we have the following corollary. 

t— >oo t 

corollary Corollary 4.1. Under the conditions of Lemma [4.41 

\nYAt) 

lim sup < 0, a.s. for each % — 1, • • • , n, 

t— >oo t 



□ 



and therefore 



In EI Yi(t) 
lim sup < 0, a.s. 

t->oo t 



Corollary 4.2. Under the conditions of Lemma [4~4 



lnfX(t)) 

lim sup < 0, a.s. for each i — 1, 

t— >oo t 



and therefore 



lim sup — < 0, a.s. 

t— >oo t 
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Proof. Recalling 



Zi(t) < Xi{t) < Yi(t),t > 0, i = 1, • • • , n 



and combining Corollary 14. 1[ we complete the proof. 



□ 



YT 



Theorem 4.4. Let the conditions of Lemma 14.41 hold. Assume further that for any t > 



and i = 1, • • • , n 



1 



Riit) := ai{t) - ^a 2 {t) + / (ln(l + 7i (t, u)) - 7i (*, u))A(du) > 0, 
and there exists constant C2 > such that 

(ln(l + li (t,u))) 2 \(du) <c 2 . 



(4.19) eql05 



(4.20) eqll3 



Then for each i = 1 , • • • , n 



(4.21) 



Proof. According to Corollary 14. 1\ it suffices to show liminf lnY ^ > o. Denote for t > 



t— >oo 



M<(t) 



d{s)dW{s) and M;(t) := / / ln(l + 7,(3, u))N(ds, du). 

Jo if 



Note that 



and by flOOD 



Since 



[M,](t) = (M i )(t)= f a 2 {s)ds<a 2 t, 
Jo 



(Mi)(t) = / / (ln(l + 7j (s,M))) 2 A(rfM)rfs < c 2 t. 
'o ./y 



1 + a 



1 



1 + s 



1 + t 



< oo, 



together with Lemma [3. 1[ we then obtain 



t— >oo t 



i-i-oo t 



JY 



lim - / ai(s)dW(s) = a.s. and lim - / / ln(l + ^(s, u))N(ds, du) = a.s. (4.22) eql21 



Moreover, it is easy to see that for any t > s 



(Ti(r)dW(r) 



and 

-t 



ln(l+7j(r, u))N(dr, du) 



o Jy 



ai(r)dW(r)- / ai{r)dW{r) 
(I Jo 



ln(l+7i(r,it))iV(dr,dtt)- / / ln(l+7i(r, u))N(dr, du). 

o Jy 
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Consequently, for any e > we can deduce that there exists constant T > such that 



ai(r)dW(r] 



< e(s + t) a.s. and 



< e(s + t) a.s. 



ln(l + 7j(r, u))N(dr, du) 

(4.23) |egl06 

whenever t > s > T. Furthermore, by Lemma [4.21 together with (14.231) . we have for t >T 



1 1 



- 1 f 

a i{ s ) - ( s ) + / (M 1 + 7i( s >«)) - li{s,u))\(du) 



ds 



+ 2e(t + T) 
+ J b u (s)exp 
+ 2e(s + t)jds, a.s. 
This further gives that for any t >T 



t . 1 /■ 

«i( r ) - ^°f( r ) + / (M 1 + Ji(r,u)) - ji(r,u))X(du) 



dr 



t . j , 

- ^i 3 ) + / (ln(l + 7i(a,u)) -7i(s,u))A(du) 

^ ./V 



t . 1 r 

a i( r ) ~ ^K r ) + / ( ln (! + 7tfau)) - 7i (?",«)) A (du) 



F,(t) " Yi(T) 

+ / 6ii(s)exp ( 
Jt v 

— 2e(t — s) — 2eTjds, a.s. 
Thus in view of ( 14.19P there exists constant K > such that for any t > T 



ds 



dr 



Hence for any t > T 



-4e(t+T) < ^ as 



- In — ^- < 4e f 1 + — ) + - In K, a.s. 
t Yj(t) ~ \ tJt 



t Yi(t) 

and the conclusion follows by letting t — > oo and the arbitrariness of e > 0. 



□ 



lemmal 



4.3 Further Properties of n— Dimensional Competitive Models 

We need the following lemma. 
Lemma 4.5. Let the conditions of Theorem 14 . 41 hold . Assume further that for i, j = 1, • • • , n 



satisfy 



iV=supij^,t>0,i^j 



Ri(t) - ^RijRjit) >0,t>0. 



(4.24) eql20 



(4.25) eqll6 
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Then 



liminf : > 0, a.s. 



t— >oo t 

where Zi(t),i = 1, • • • , n are solutions of (12. 5p . 



(4.26) 



Remark 4.1. For i,j = 1, 
condition (I4.24p must hold. 



. n 



and t > 0, if %(£) takes finite-number values, then 



Proof. It is sufficient to show lim sup \ In -jK-r < 0. Note from Lemma 14.21 that for any 

t > S > 
1 1 



exp 



Zi(t) Zi(s 

o-i(s)dW(s 
+ / &n(r)exp 



fl i( r ) - y]M r ) y }( r ) - o a t( r ) + / Qn-0-+li(r,u))-'yi(r,u))\(du) 

■ , ■ A JY 

1^3 



dr 




ln(l + 7j(s, u))N(ds, du) 



a i (T)-J2b ij (T)Y j (r)--af(T) 

i¥=3 



+ / (ln(l + 7i(r,w)) -7i(r, u))X(du)^dr 



a^dWir) 




ln(l + 7j(r, u))N(dr 1 du) )dr. 



(4.27) eqll5 



Applying the Ito formula, for any t > s > 
J bu(r)Yi(r)dr = \nY t (s) - InY^t) 
+ 



ds (4.28) eql07 



a i( r ) ~ o°f( r ) + / (ln(l + 7<(r,«)) -7i(r,«))A(tZ«) 
^ </y 

+ / <7i(r)<W(r)+ / / ln(l + 7i(r,u))iV(dr,du). 

Js </s </Y 

This, together with Theorem 14.41 and (I4.23p . yields that for any e > there exists T > 
such that 



bii(r)Yi(r)dr < 



a i( r ) ~ 2°K r ) + / ( ln ( 1 + u)) - 7i(r, w))A(tZ«) 



(4.29) eqll4 



+ 3e(s + t) 
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whenever t > s > T. Moreover taking into account (14.281) and (14.291) . we have for t > s > T 



bij(r)Yj(r)dr 



bij(r) 



ds. 



b n( r ) 

< Rij [ hAr)YAr)dr 

J s 

< 3e(s + t)Rij 
+ J 

Putting this into (I4.27P leads to 

Yjt) = z^y exp ( ~ / N r) ~ E R ^ R ^] dr + < s + *) ( 3 E ^ + 2 

[^(t)-^^-(t 
+ e(r + t) (3j2 R ij + 2 )) dr , 
which, in addition to (I4.25p . implies 

Y^t) = Y{s) exp ( e(s + 1) ( 3 S ^ + 2 )) + [ bli (r) exp ( e(r + 1} ( 3 S ^ + 2 ) ) rfr - 



+ / 6jj(r)exp 



Carrying out similar arguments to Theorem 14. 4[ we can deduce that there exists K > such 
that for t > s > T 



exp ( - 2e( S + t) (3 £ fly + 2) ) -L- < K 



and the conclusion follows. □ 

Now a combination of Theorem 14.41 and Lemma 14.51 gives the following theorem. 

Theorem 4.5. Under the conditions of Lemma 14.51 for each i — 1, • • • ,n 

, lnXJt) 

lim = 0, a.s. 

t— >oo t 

Another important property of a population dynamics is the extinction which means 
every species will become extinct. The most natural analogue for the stochastic population 
dynamics ( 11.51) is that every species will become extinct with probability 1. To be precise, 
let us give the definition. 

Definition 4.2. Stochastic population dynamics (I1.5P is said to be extinct with probability 
1 if, for every initial data Xq G M™, the solution Xi(t),t > 0, has the property 



lim Xi(t) — > a.s. 

t— >oo 
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Theorem 4.6. Let assumption (A) and (I4.20p hold. Assume further that 

1 /■* 

7]i := limsup - / (3i(s)ds < 0, 

t— S-oo t Jq 

where, for t > and i — 1, • • • , n, 

:= Oi(t) - ^o-f(t) - f ( 7f (t, w) - ln(l + Ti (t, w)))A(dw). 



2 

Then stochastic population dynamics ( 11. 5ft is extinct a.s. 

Proof. Recalling by the comparison theorem that, for any t > and « = 1, • • • , n, 

Xi(t) < Yi(t), 

we only need to verify limsup^^ Yi(t) = a.s., due to 

< liminf X t (t) < limsupX^t) < limsup 

Since bi(t) > 0, by (14 .7p it is easy to deserve that 

Yi(t) < Xi(0) exp ( f Pi{s)ds+ [ o-i(s)dW(s)+ [ [ ln(l + 7< (s, u))N(ds, du)) 
^ Jo Jo Jo Jy ' 

= X,(0)exp (t(]J Pi(8)d8 + jj <Ti( 8 )dW(3) 

+ jjf ln(l+ 7 i(s,w))iV(ds,^))y 

Thanks to rji < 0, in addition to (I4.22p . we deduce that lim sup^^ Yi(t) = a.s. and the 
conclusion follows. 

□ 

Remark 4.2. In Theorem 14.31 we know that one dimensional our model has a unique 
invariant measure under some conditions, however we can not obtain the same result for 
ndimensional model (n > 2). 



5 Conclusions and Further Remarks 

In this paper, we discuss competitive Lotka-Volterra population dynamics with jumps. 
We show that the model admits a unique global positive solution, investigate uniformly finite 
p-th moment with p > 0, stochastic ultimate boundedness, invariant measure and long-term 
behaviors of solutions. Moreover, using a variation-of-constants formula for a class of SDEs 
with jumps, we provide explicit solution for the model, investigate precisely the sample 
Lyapunov exponent for each component and the extinction of our n-dimensional model. 
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As we mentioned in the introduction section, random perturbations of interspecific or 
intraspecific interactions by white noise is one of ways to perturb population dynamics. In 
[13], Mao, et al. investigate stochastic n-dimensional Lotka-Volterra systems 



dX(t) = diag(Xi(t), • • • , X n (t)) [(a + BX(t))dt + aX(t)dW(t)} , (5.1) eq53 



where a = (a 1; ■ •• ,a n ) T ,B = (bij) nxn ,o~ = (c7y) nxn . It is interesting to know what would 
happen if stochastic Lotka-Volterra systems fl5.ll) are further perturbed by jump diffusions, 
namely 

dX(t) = diag(Xi(r), • • • , X n (r)) \(a + BX(t))dt + aX(t)dW(t) 

(5.2) I eq54 



+ / 7(X(t~),u)N(dt,du) 
Jy 

where 7 = (71, ■ • • ,7 n ) T - On the other hand, the hybrid systems driven by continuous-time 
Markov chains have been used to model many practical systems where they may experience 
abrupt changes in their structure and parameters caused by phenomena such as environ- 
mental disturbances [15]. As mentioned in Zhu and Yin [251 [26], interspecific or intraspecific 
interactions are often subject to environmental noise, and the qualitative changes cannot be 
described by the traditional (deterministic or stochastic) Lotka-Volterra models. For exam- 
ple, interspecific or intraspecific interactions often vary according to the changes in nutrition 
and/or food resources. We use the continuous-time Markov chain r(t) with a finite state 
space M. = {1, • • • , m} to model these abrupt changes, and need to deal with stochastic 
hybrid population dynamics with jumps 

dX(t) = diag(X x (r),--- ,X n (t-))\(a(r(t)) + B(r(t))X(t))dt + a(r(t))X(t)dW(t) 

(5.3) I eq55 

+ / j(X(t~),r(t),u)N{dt,du)' 



We will report these in our following papers. 
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